clear,clc;
x=fsolve('fun_exp1',[10,10]',optimset('Display','off'))
